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Abstract 

Spontaneous symmetry breaking in (l-l-l)-dimensional 0* theory is studied with 
discretized light-front quantization. Taking effects of non-diagonal interactions 
into account, the first few terms of the commutation relations [oo, an] are recal- 
culated in the H expansion. Our result of the critical coupling is still consistent 
with the equal-time result 22^^/?i < Acr < 55.5^^/?i. We also have examined 
effects of regarding the ratio of the bare coupling constant to a renormalized 
mass as an independent parameter in the h expansion. 



1 Introduction 

Light-front QCD appears to be a hope of understanding hadrons from a field- 
theoretical point of view.^) Light-front field theory (LFT) is expected to be equiv- 
alent to the ordinary field theory not only in the perturbative region, ^-^ but also 
in the non-perturbative region. Spontaneous symmetry breaking(SSB) is a typ- 
ical phenomenon in the non-perturbative region. One of the most important 
feature of LFT is that the true vacuum is identical to the Fock vacuum. There- 
fore it appears that SSB does not occur in LFT, but it is now widely believed 
that a rich low-energy structure of the theory can be explained through zero 
modes of field operators. In LFT, zero modes of some field operators are not 
independent variables, but obey certain constraint equations and are given by 
complex functions of non-zero modes. In the ordinary equal-time formulation, 
(l-l-l)-dimensional (/>^ model exhibits spontaneous breakdown of Z2 symmetry.^) 
It is therefore an important issue to confirm that the phase transition of (j)'^ model 
is treated correctly in the light-front formulation. Several authors^' have in- 
vestigated the phase transition of the cj)'^ model in the light-front formulation. 
In this formulation a non-zero vacuum expectation value of the zero mode is a 
signal of the phase transition. In the previous paper^^ the author and Yahiro 
have pointed out that if we can calculate the commutation relations of the zero 
mode and the non-zero modes [ao,a„], we can compute the critical coupling of 
the (1 -I- l)-dimensional light-front (j)* model. In Ref. 9) we calculated the main 
part of [ao,a„] to 0{h'^) in the h expansion. We neglected non-diagonal inter- 
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actions, which cause transitions among different modes, since they only give 
vanishing contributions for certain vacuum expectation values at the first level 
of the interactions. In this paper wc take the non-diagonal interactions into 
account besides diagonal interactions, since the non-diagonal interactions give 
non-vanishing contributions at the second level of the interactions and subse- 
quently the critical coupling will shift. Taking all the interactions into account 
we evaluate the critical coupling within the third-order corrections of the h ex- 
pansion. We also calculate a renormalized mass using front-form pertubation 
theory and compute the critical coupling in terms of the renormalized mass. 

This paper is organized as follows; In § 2 we denote some conventions of 
light-front theory and formulate (1 -|- l)-dimensional cj)'^ model. In § 3 we calcu- 
late the first few terms of the commutation relations [oo, ffl,i] in the h expansion. 
In § 4 we calculate the vacuum expectation value of the zero mode and compute 
the critical coupling. In § 5 we calculate the relation between a bare mass and 
a renormalized mass and evaluate the critical coupling in terms of the renor- 
malized mass. In the last section we give summary. The text is supplemented 
by five appendices. In these appendices, several relevant commutation relations 
are given to calculate [ao5an]- These commutation relations are related each 
other. 



2 The model and formulation 

We define the light-front coordinates a;^ = (a;° ± x^)/V^. The Lagrangian 
density of (1 + l)-dimensional (f)'^ theory is written 

£ = 5+</.5_<^-^<^2-A<^l (1) 

We put the quantum system in a box of length d and impose periodic boundary 
conditions. The field operator (f) is then expanded as 

cl>{x) = -^ao + J2 -j^an{x+)e<^- , (2) 
v47r ^„ y^4TT\n\ 

where A;+ = 2^12/(1. The coefficients of the expansion are operators which satisfy 
the canonical commutation relations, ^^'''^'^^ 

[ak,ai] = [al,a}] ^ 0, [a*;, a|] = Mfe,;, k,l>0, (3) 

where al. = o_fe. The zero-mode operator oq is not an independent quantity, 
since it governs the constraint equation^) '''^ 



-|-a„aoajj + alaoan - 37iao) + 6X^3 = $(a, a^), (4) 
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where g = 247r/i^/A is the only parameter of the theory and we have added the 
term J2n>i "TT^^o *° remove tadpole divergences and J^n defined as 

The commutation relations and the constraint equation were obtained 
by quantizing the classical system with the Dirac-Bergmann quantization pro- 
cedure. This procedure does not specify an operator ordering in the constraint 
equation (Q). Hence we have assumed the Weyl (symmetric) ordering. -'^^^ 

In the previous paper^^ the author and Yahiro have emphasized the im- 
portance of the commutation relations [ao,a„] in calculating the vacuum ex- 
pectation value (0|ao|0). We obtained the first few terms of [ao,a„] in the fi 
expansion, neglecting the cubic term in (^) which causes the non-diagonal 
interactions for simplicity. In fact, at the first level of the interactions, it does 
not contribute to the vacuum expectation value (0|$(a, a'l')|0). We also used 
[ao,a'^am] = 0, which hold when the cubic term is absent. In this paper we 
recalculate (0|$(a, a'l')|0) within 0(h^) without neglecting the cubic term. The 
cubic term gives non-zero contributions at the second level of the interactions. 



3 Commutation Relations [ag, a^] 

In this section we will calculate [uq, a„] from 

[$(a,at),a„] -0. (6) 
Using the commutation relations (^), Eq.(^ is written explicitly as 

A[ao, an]\-aoanh + -[ao, a„]?i|- — [[[oo, a„], a™], aj„] 
n n b ^ — ' m 

m>l 

+ ^ao[[ao,an],ao] + ^[[[ao, a„], ao], ao] = 0, (7) 



where 



m 

m>l 



In this paper, we will use the letter m as a natural number and the letters k and 
I as integers. We mainly use the letter n as a natural number except for some 
cases, where we use the letter ti as a non-zero integer. The vacuum expectation 

^The operator A in Ref. 9) is diflferent from A (1) by a factor 6. 
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value of the fifth term(~ ^ akcin-k) in (0) which originates from the cubic term 
^3 in (Q) vanishes and the last three terms in vanish or become higher order 
in ao as long as the fifth term is neglected (oq is very small around the critical 
point ). In this reason in the previous paper^^ we have neglected the last four 
terms in (^. 

We will fix [flg, a„] order by order in the h expansion. At the lowest level we 
have 

[ao,a„] = ^aoa„S,+ -5- ,,, ^ =afcan^fc?t + 0(?l^). (9) 

The right hand side of the commutation relations [oq, a„] can be divided into 
two groups by the powers of the zero mode oq. We call such terms as the first 
term in diagonal that contain odd powers of ao. In the contrast, we call such 
terms as the second term in non-diagonal that contain even powers of oq. 
Note that the powers of oq can increase only through the commutation relations 

Since there is no a priori way of fixing the operator ordering in (p|) , we have 
fixed the operator ordering so that we can calculate the vacuum expectation 
value easily. Although the difference of the operator ordering affects the higher- 
order terms in the expansion, the vacuum expectation values (0 Jao, ao]ajj|0) 
would not be affected. We get [ao,aJJ by Hermite conjugation of (1) 

[ao,at] = -aiao^n - ^£ ^^^^^^ "^^"^ 2^^ + ^^^'^ 

(10) 

Owing to the difference of the operator ordering, the 0{h^) term in the second 
line differs from that in the first line. We adopt the operator ordering in the 
second line, where A and oq are rearranged in order from the left. By noting 

at = 



a_„ and changing some variables as k ~> — fc, (^) and (|To|) can be combined 



as 



11 1 

[ao,a„] = — -aoa„?i+ — ^=====afca„_fe?i + 0(?i^), n ^ 0, (11) 



where e(n) — l{n > 0),e(n) — — l(n < 0). Note that the second term in ( |ll[ ) 
are quadratic in the non-zero modes due to the non-diagonal interactions. 

Before proceeding to higher order in h, let us examine what terms should be 
left in our approximation. We consider the following expansion (see Ref. 9)) 

[ao,a„] EE ^[ao, a„]p?i^ 
p>i 
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p>l p>3l>l 

+ (higher — order terms) + ©(oq), n > 1, (12) 

where ap{n), (3p^i{n) and 7p,fc; are the coefficients of the diagonal terms to be 
determined. Now we comment on the the expansion (|l2|). Since a;]9 -^aih, 
[ao,a; 3 ^aoaifi and the same relations hold for aj, the terms -^ajaiaoanfi^ 
in (|l2) are produced from the fourth term in (^. The non-diagonal terms do 
not directly give non-zero contributions, but they affect the diagonal terms at 
higher orders, and at 0{h^) we will write down [ao, a„] completely. At 0(?i'^) and 
0{h'^) only the terms aoanhFip — 3,4) give non-zero contributions to the 
vacuum expectation value (0|<i>(a, a'l')|0) within 0(?i^). For example, the terms 
^^^4"^ a\aiaQanh^ can indirectly give non-zero contributions to (0|$(a, a^)|0) 
only at 0{h^). At 0{h^) and 0{h'^) therefore we only pay attention to terms 
that reduce to aodnh^ or aoa„?i^. 

At 0{h^) only four terms in the constraint equation (Q) are concerned. These 
terms are calculated from [oq, a„]i and a„]i. The fundamental commutation 
relations [^,a„]i and some secondary commutation relations are denoted in 
Appendices A, B and C. From (^ij), (|5|) and ( |6l|) we obtain 

[ao, On] 2 

1 1 
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ao{aian-k-iak + akaian-k^i) 
^^\kni^t:-k {n-kW\nkl{n-k-l)\ 

4 ^ Vn k n-kJ\A^ A^ 



k^O.n 



k)\A^ A^ V ^\nk{n-k)\ 



za^an 



e(m) / I ^ "A \^ o.n{a-makam-k — ima-ka-m+k) 
Z^^^Ka^- II N ^2^^ ^\mk{m-k)\ 

^ ^ ^ ^'■"^^ 



m=i/^o,nfe#o,m ' l)km{m - k)\ 

X ao[a„_;ai, a_mafcam_fc - ama_fcafe_m], 727^0. (13) 

Although the last term in ([l3|) belongs to [ao,a„]3, we write it down here so 
that we may forget it at 0{h ). We have left the symbols e(m)(= 1) as they are, 
since they are useful to distinguish the non-diagonal terms from the diagonal 
terms. 
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At 0{h^) six terms in the constraint equation are concerned. To obtain 
[ao,a„]3 several commutation relations at 0(h) and 0{h^) are needed. Rele- 
vant commutation relations are denoted in Appendices. Non-diagonal terms 
in [flo, a„]3 do not give non-zero contributions to the vacuum expectation value 
(0|$(a, at)|0) within 0{h^). Neglecting irrelevant terms, from (|l3|) , (|5|) , , @ 
and ( |67| ) we have 

( 4 ^C(2)\ 1 
1 ^ ^ / 1 11 5 



^4 ^ ^ V6n 48fc 48Z 8(n - /c) 48(n - A: - 

1 , , 

ao(aian-k-iak + akaian-k-i) 



(n — k) y''\nkl{n — fc — /)| 



(44 X! X! 



^774 X! 2 / \ ( \ui =aoan-,nakam-k 

m=l(m^n) fe#0,m V I v / \ /I 

^7944 X! X! 2 / i; ^ :aoa;a^_;an-^ 

-L^^ 1 / ,77r TO^v — ^ "11 



12^' .^^n.^fc-. (" - kW\knlin-k-l)\ 

- + -j^akaian-k-i^, n > 0. (14) 

We only have to pick up terms including aoa„?i^ in (Q) to calculate [oq, a„]4. 
From (H) and (|6^) we have 

/ 7 5C(2) C(3)\ 1 
[ao, a„]4 = - + ^ + j ^aoa„. (15) 

4 The critical coupling 

Using the commutation relations [ao,a„] in the previous section, we calculate 
the vacuum expectation values (0|[ao, a„]4ajj|0) to compute the critical coupling. 
Sandwiching (||) with (0| and aJjO), we have 

{O\[ao,an]ial\0)h^^-ah^, (16) 
q n 
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where a = (0|ao|0) and g — g + 3cr^. In the same way from (|l|) we have 

'6x4/ 7 5^(2) ((3)^ 



[ao, a„]4ajj|0)?i'' 



gJ \3n^ Qin? 



2n 



ah" 



(17) 



where C,{n) is the Riemann zeta function. Using identities such as X]fc=i k(n-k) " 



I Till from @ and (0) we have 

'6\2 1 



(0|[ao,a„]24|0)?i' 



()■ 



ah'' 



" fc=i 

n-l 



.g/ V2n2 ^ fc2 + 2n3 ^ k)'^ ' 



(18) 



6\4/J ^ iLvl ^LV— h^ 



k—l k—l 



(19) 



The 0{h'^) and 0(?i^) terms in (|l|) and O(n^) terms in (|l|) originate from the 
non-diagonal interactions. CoUecting the above resuhs, we have 

{0\[ao,an]al\Q) ^^-ah^ - \ah^ 
gn \g/ 



6\3/ 4 



C(2) 



1 



^ fc=i 



1 "^^ 1 

„2 h ) ''^^ 



6\4/ 7 



() 



5C(2) , C(3) 



1 



C(2) 



5/ V3n4 6ri2 2n 3n4 



7 
3^ 



ri-l 



^1 ^Y-J_\ .5 

fc + 12„2 pj'^'^ • 

fe=l fc=l 

Sandwiching (||) with the Fock vacuum |0), we have 
00 ^ 

=a^+ga-y2-{0\[ao,a,,]ai\0) 



(20) 



n=l 



^ ^ ■-' ri=l fc=l 
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(•AO n oo _ n— 1 oo _ n-1 



n=l fc=l n=l fc=l 



(21) 



At the critical point non-trivial solutions start to appear besides the trivial 
solution (7 = 0. From ( ]2l| ) the critical point is given by the condition that the 
linear terms in a balance to zero. ©This condition is written as 

= 36 - a;2(l - 0.2326a; + 0.1665a;2 - 0.1068x3), (22) 

where x = Gnh/g = and we have substituted the values X^J^i ^ X]fc=i i 

0.3529,Er=i_Y w = 0-2288, 

TIT X]fe=i W ~ 0.0966. Owing to the non-diagonal interactions, the third 
and the fourth terms in the parentheses in ( p^ have increased by some 15% 
and 8% respectively over the previous results in Ref. 9), where the non-diagonal 
interactions are neglected. 

We adopt the Pade approximations^^' to the alternating series 1 — cix -\- 
cix^ — c-iX^ + • • •. This series is approximated as [Af/A^], which denotes a 
rational equation of polynomials of degrees M and N: 

l-cix + --- + (-l)'e. , (23) 

l + (^^-ci)x 



1/21 . (24) 



Our results are 



[0/0] Aer = 2V/?i, [1/1] A„ = 28.lAiV;i, 

[0/1] A„ = 46.0AtV7i, [2/1] A„ = 26.0ii'^/h, 

[1/2] A„ = 25.5mV^- (25) 

In Ref. 9) we have assumed that the alternating series has such properties that 
the lower (upper) bounds of Acr are obtained from cases M + A^=even(odd). 
The values in (|2q ) exhibit that this assumption does not hold. However, all 
the values in ( |25| ) still lie within the range 22fi^/h < Acr < 55.5/i^/?i of the 
equal-time result. 

The ratio of the bare coupling constant to the bare mass is the only pa- 
rameter of the theory and in the literature gcr is computed. In two dimension 
the coupling constant and the mass suffer only finite quantum corrections in 
this model. However, in four dimension the bare quantities become infinite and 
it will be adequate to describe the theory by renormalized quantities. Hence 
as a preliminary to study four-dimensional models, we adopt the ratio of the 
bare coupling constant to a renormalized mass as a new parameter in the next 
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section. We avoid treating coupling constant renormalization since it seems not 
so easy to calculate quantum corrections to sufficient order in Ti. In the next 
section we compute the critical value of the bare coupling constant in terms of 
the renormalized mass. 



5 The critical coupling in terms of the renor- 
malized mass 

Let us study a relation between the bare mass fi and the renormalized mass m 
for the model. We split the Hamiltonian into two parts H = ^^{Hq + Hj), 

where Hq = 5X^2 ~ ffS^i and Hj = GJ^a^^zm- The Hamiltonian 
HzM contains all of the zero-mode interactions 

TT ^0 \ ^ / 2 \ 

nzM — — T + T/ 1— rlanflo'^-" ~ '^oOna-naoj 
4 4 ^-^ \n\ 

{auaoaian + akaia^^a^). (26) 



Hal\0) « ^^hal\0). (27) 



The states |0) are eigenstates of the unperturbed Hamiltonian Hq with eigen- 
values = ^^Y~?i.@We approximately define the renormalized mass m through 
eigenvalues of H as 

247rm2 , 

If the interactions were absent we would obtain the bare mass /i, hence this 
definition of the renormalized mass m will be plausible. We will solve the 
eigenvalue equation ( p7|) perturbatively. 

Effects of the interaction Hj have already been calculated to the second 
order in perturbation theory in the limit c? — > cx) in Ref. 8). Taking ( p2[ ) into 
consideration, we calculate effects of the interaction Hi to the third order in 
perturbation theory in the limit d — > 00 . 

First, we calculate the contributions of the non-zero modes. Since we have 
added the counter term in (^) , the tad pole diagrams always vanish. At the sec- 
ond order in perturbation theory we have a non- vanishing diagram (see Fig. 1): 

(0|«p(6E4) p_^^ (6E4)4|0) 

i+j<p 

_ A ^ 



9P ^j^o iPiP-^-3){^ 



( dx ( dv- -Ti^ — Ti^ , 

gpJo Jo ' {x + y){P-x){P-y) gp A 



(28) 
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where i,j and p are natural numbers and E = is the unperturbed energy; 
we have set x = 2m/d,y — l-Kjjd and P = 2'!Tp/d and have taken the Hmit 
d, p ^ oo so that P remams finite. Equation ( [2^ is independent of d and 
it holds also at large d. Therefore in the following we consider the situation 
d ^ oo. In the ordinary equal-time formulation the corresponding Feynman 
diagram of Fig. 1 gives the same result. In the same way at the third order in 
perturbation theory(see Fig. 2) we have 

(o|ap(6E4) p ^^ (eEJ p ^^ (6E4)4lo) = (29) 

where / is the following quantity which approaches a definite value in the limit 
d ^ oo 

1 = 

k<i+j<p ^ 

^ iv(v i ^ ^ ^ ]k(i \ j k)( ^ ^ ^ ^ ) 

i,j,k>0 '^^^ J'\p i j p-i-j I y 'J '\p k i+j-k p-i-j / 



P nP-x .x+y P'^(P~X~y) 

dx dy dz- -^-r— — — — = 1.645. 

' ^ ' {x + y)^{P - x){P - y){P - X - y + z) 



Jo Jo 



(30) 





P-i-J P-i-J 



Fig. 1. The second order diagram. Fig. 2. The third order diagram. 



Second, we calculate the contributions of H^m- Solving the constraint equation 
(||) perturbatively in 1/g for ao, we have 

ao = --Y.3 + "T X!(S3»n"-" + ana_„E3 + ""E30-" ^ IEs) ^ ■ (31) 
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Substituting this expression into Hzm (p6[), to 0{h ) we find 



1 1 



9P Hp -k) ^ k{p + k) 



^ 18 1 ^ 



where we have neglected some constant terms, which means a shift of the vacuum 
energy. Note that X^Ci k(p^-k) ^^^"^ I]fc>o fc(p+fc) ^^'^ 0{Ti^) terms vanish in 
the hmit p oo. The 0{h^) terms also vanish in this hmit. Thus the mass will 
not be corrected by the zero mode in this limit. However, the convergences in 
( ^8| ) and ( |30|) will be improved by inclusions of the zero mode.*-* 

The eigenvalues of Hq and their corrections ( |2^ ) and (|2^) have the common 
factor 1/p as expected, and we get the following relation between the renormal- 
ized mass m and the bare mass 

14.80^2 88.83. 
9 9 

Now we redefine the parameter g as g — 247rm^/A, which we regard as indepen- 
dent of h, instead of g = 247r/x^/A. In the following g means this new parameter. 
From (p3|) we have 



/i'(i —h^ + —^n'^) + o{n^). (33) 



247r^^ 14.80^2 88.83^3 

— 7 — =.9+ ^ 2 

A 9 9 



'^h^ + o{h'^). (34) 



Taking this Ti dependence into account, the constraint equation (|^) becomes 



l.h.s. of (7) + ^^'^^ [ao,an]h^ r^[ao,a„]?i^ = 0. (35) 

Off Off 



14.80. , o 88.83 
— - — [ao,a„J?i 

Off 

These two terms affect [ao,a„]2 and [ao,a„]3 

[oo, a„]3 = r.h.s. of (7) - ^^'^!^ -Qoa„, (36) 

6ffy4^ n 

^ 14.80 1 88.83 1 

[ao, a„j4 = r.h.s. of (7) + ^^^^J^^oOn + ^^2^42 ~^o""- (^'^ 

Thus the fluctuations in the ?i expansion of [ag, a„] have diminished in terms of 
the renormalized mass. This time we have instead of ( p^ ) 

= 36 - x^{l - 0.23261X + 0.124889a;2 - 0.0742a;3), (38) 
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with x = ^. 

Using again the Pade approximations, this time we have 

[0/0] Acr==24TOV^, [1/1] A„ = 29.7mV?i, 
[0/1] A„. = 46.0mVn, [2/1] Acr = 33.3mV?i, 
[1/2] A„ = 32.2mV?i- (39) 



The variance of the values in ( |39D is smaUer than that of the values in (|2^) . The 
values in ( ]39| ) seem to approach their mean value as the order Af + of the 
Pade approximations increases unlike the case 



6 Summary 

Taking the non-diagonal interactions into account besides the diagonal interac- 
tions, we have calculated the relevant terms of [ao,a„] to 0(?i^). In addition 
to the Riemann zeta function ^(s) = 'Y^=\ ll^^ more complicated fmiction 
SJ^i W appeared corresponding to the non-diagonal interactions. 

Owing to the non-diagonal interactions, the third and fourth coefficients of the 
alternating series in (^2|) have increased by some 15% and 8%. We have applied 
the Pade approximations to the alternating series to compute the critical cou- 
pling. Although the critical coupling shifts sensitively to the coefficients in the 
h, expansion, at each order of the Pade approximations Acr hes within the range 
Tl[i} jTi < Acr < 55.5fj,^/h of the equal-time result. We also have calculated the 
renormalized mass in front- form perturbation theory. We have obtained another 
h expansion of [ao,a„] by regarding the ratio of the bare coupling constant to 
the renormalized mass as a new independent parameter. We have found that 
the fluctuations in the h expansion of [ag, a„] and the variance of the values of 
the critical coupling have diminished. 
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Appendix A( 



A' 



a. 



In this appendix we calculate the fundamental commutation relations [^,a„]. 
We can obtain [;j,a„i] from [00,0^] as follows 



If .1 1 



12 



( 



,4\ 



n \ 1 

m J A 



m \m/ A-^ \niJ 

73 X! 



1 



/ m. 

2 



2 * 
^3 ^Q^m'^ 



2A3 



■\/|fcm(TO — fc)| 



akajn^kaoh + 0{h ), m > 0. 



(40) 



The abbreviated terms lead to higher-order contributions. In the same way we 
have 



r 1 

J, a. 



h 1 

771 



— y 



?i\2 1 

TTlJ A^ 

1 



?i\3 1 



1 



2 t R 
^3 '^O^m" 



^|fcm(TO — fc)| 



aL-fcaI«0?i + O(?i^),TO > 0. 



(41) 



Equations (|o|) and (|l|) are combined as 



A' 



h 1 /^y i /ny i 

n Vn/ ^3 + V„J 



'?i\2 1 



n\3 1 



nA3 



e(n) 1 



akan^kO-a^ + 0{h'^), n ^ 0. 



(42) 



Appendix B 



In this appendix we calculate the first few terms of [ao, al^am] in the h expansion. 
In our approximation it is sufficient to get [aQ,al-,^am]p{p — 1,2,3).@ More 
precisely, we need the coefficients of the terms in the expansion 



[ao, al^am] — {aoaah + aoaah^ + aoaaaah^ + aoaah'^) + (aaah), 



(43) 



where the parentheses denote the disregard of the coefficients and a's represent 
the non-zero modes. We have neglected terms such as (apaaTi^) and iaaah^) 
that become higher order in fi. 

Let us examine how [ao,aj„am] would be if the non-diagonal interactions 
were absent. At the lowest level we have 



[ao,aLa 



mjl 



ala[ao-,a.m]i + [ao,a]',Jia„ 



A 



aln [a''iai,ao]a„ih. 



(44) 
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This equation will indicate that if we neglect the non-diagonal interactions, 
[ao,al^am] can vanish at any order of h. In the following we calculate the 
coefficients of (|4^) taking the non-diagonal interactions into account. 

First, we calculate the contributions of the diagonal parts [ao, am\dp 
(p = 1, 2, 3). Using (|n]) and (|l|), we have 

2 3 

J\ in Ji Tn 

® {aQaah^ + aaah^ -j- a^aaali? + aoaaaah'^). (45) 

Here the symbol ® means that the following terms can be neglected in our 
approximation. In the following we usually neglect such higher-order terms 
that become 0{h^) in (0| [ao, a„]ajj0). The 0{h) terms have canceled, since 
— — -I- — = 0. In the same way we have 

mm ^ 

[ao, aL]d2am?i^ + aJ„[ao, am]d2?i^ 

^ t fc2 3 t t;3 

— jaoa'amTL 7^5 — ^-aoa' am?! 



1 

ao{al^aiam-k-iak 



1 1 



fc^CmW-fe ("^ - k)^\mkl{m-l-k)\ 

+ al^akaiam-k-i + alal^_f^_ialara + alajal^_f^_iara)fi^. (46) 

The last term derives from the first e^-term in (^3|). At 0{h^) we only have to 
pay attention to the terms including aoaJ„am?i^. Since— + = = 0, 

from the first term of (|lj) we have 

[ao, aln\d3amfi^ + aJ„[ao, am]d3?i^ = 0?i^ ® (aoaooo?i^). (47) 

Second, we calculate the contributions of the non-diagonal parts [oq, o,l^am]ndp{p - 
1,2). At the lowest level we have 

[ao, a,J,Jndiam?i + ajjjao 

— 0,1,, — 7 — ajfi jh 

— 7 / — , a^ha—m+kO'miT' 

2^ M0,m Vl^"^("^-fc)l 

e2(m) V- 1 .2 



2 Vm^ 



/\u^ i\i ^ =aoa-;a;_TOafcaTO~fc^^ 



(48) 
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Since [^aaah^,an] = {aahF^^) + (aoaaaaah^^^), the contributions of the first 
and third terms in (Esh to (0| [oq, a„]ajj|0) are 0{h^) and 0{h^) respectively. 
The second term in ( [48[ ) gives non-zero contributions to [ao, amaj^jds- Hence 
we cannot neglect the second term. In the same way we have 

[ao, a],Jnd2am^^ + aJ„[ao, am]nd2^'^ = ®{aaah^). (49) 

The right hand side of ( |49| ) gives only higher-order contributions. Thus as 
concerned with the non-diagonal interactions, we only have to take the lowest- 
order terms into account . 

Collecting the above results for the diagonal and the non-diagonal interac- 
tions, we obtain the main result in this appendix: 

_}_ V- \ 

i^h-k (- - kWl^klim -l-k)\ 
X aoial^aiajn-k-iCLk + al^akaiava-k-i 

^ ^ =aoa_;a_„i+;afcam-fc?i^.@ 



kk^n ikr^. m^mn.~l){m-k)\ 



(50) 



Using (50), we can write down the lowest term of the expansion of [[ao, a„jajj^], a„] 
as 



[[ao,a™a]^Ji,a„]i?i^ 
_e(n) V- \p _ 

4yl3 



e(m) 



i^Q^nk^o,m yj\ln{n~l)km{m-k)\ 

X aaan^iai{al-fika7yi-k - o.malal^_^)h^ 

e(n) / 1 1 1 \ 1 1 



E 

k^O.n 



n k n-kJ y/\nk{n-k)\ 
e(m) 

fc^o.m V\^Hm-k)\ 



afca„_fc?i 



1/1 1 2\ ^l*^) ^ t t t ^■^2 



(51) 



The first term becomes 0{h^) when sandwiched between (0| and aJjO). The 
second term leads to 0{h^) contributions through ao[[ao, a„], ao](see (60)). The 
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next higher-order term [[ao, amaj„]i, a„]2?i'^ produces only higher-order terms. 
The commutation relations [[ao, al^am]2, an]i contribute to [ao, anjs in the form 
(aoaaah^): 



— [[ao,a™a|„]2,a„]i?i' 
m 

m— 1 



{n - A;)v/|nfcZ(n-/c-/)| 



xao{aia„-k-iak + akaian^k-i)'h^ 



' /^"^^ mfm — k)\/\mkn(m — fc — 7i)| 



A3 



2A3 



m=l(m5iin) k^O,m 
4^3 E E 



m2-\/|n(m — ri)k{m — fc)| 



1 



O'OO'n-mO'kO'm-k^ 



(^oa'ka^_n-k<^nJ 



/ s, ,n m(m — n)\/\mkn{m — k — n), 

m,—l[m^n} k^O,rn — n ^ /VI V /I 

® {aoalj^aa + aoaaam+n + ■ ■ ■)fi^ ■ (52) 

The terms in the last two lines become higher order, since amoJilO) = 0(to ^ n) 
and (0|aj„ = am+„ajj|0) = 0, etc.. 



Appendix C ao[[<^05 CLn] , ^o] 

In this appendix we calculate ao[[ao, a„], oq]. Expanding [ao,a„] in h, we have 

ao[[ao, a„], ao] = ao[[ao, a„]i, ao]?i + ao[[ao, a„]2, ao]?i^. (53) 

We calculate these terms by dividing [ao,a„]p(p — 1,2) into the diagonal and 
non-diagonal parts. 

First, we calculate the first term in (p3[). The contribution of the diagonal 
part of [ao,a„]i is calculated as 

ao[[ao, a„]di, ao]?i = ao 



1 

ao- 

n 



nA 
1 



aofln, ao 



aoa„?i - ao^ao[ao, a„]?i 



®(aoaaaaoa?i^) ® ^ao-^aoa?i^^ © @^ao-^aoaa?i' 



e(aoaaaaoa?i^) G ^ ^ aoa„?t^ 



'2nA^ E 



e(n) 



„ ^J\kn{n-k)\ 



(54) 
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where we have used [^,ao] ~ aaah. The three terras in the last hne only 
produce higher-order terms. Thus within our approximation the contribution of 
the diagonal part of [qq, an] is negligible. @The contribution of the non-diagonal 
part of [ao,a„]i is calculated as 



ao[[ao 



A 



ako-n-k — -^afc[ao, a„_fe] 



(55) 



The first term in the parentheses in ( |55| ) is calculated as 
^ <n) r 

t{n) NT- NT^ 1 



AA^^\kn(n~k)\ -VIM- "01 



-amoj'a|„_;)afea„_fe?i^@ ® (aoaaaaa?i^). (56) 
At the lowest level the second term in the parentheses in (p5|) is calculated as 



6(n) 



■y/|/cn(n — fc)| 



-ao 



— y 



e(n) 



-.^n (" - fc)v/|M"-fc)r''^ 



1 fi2 

0-k—raoan-kn 



-ao 



- y 

OA ^ 



2 A .^^^ ^\kn{n~k)\ 



E - 

l=^0,n-k 

-y 

A2 



e(n ~ /c) 



1,, 2A ^|;(n-fc)(n-fc-0| 



aia^-k-fi 



e{n) 



2^2 



2 t:2 

aQakan-kii 



4^2 E E 



e(n) 



{n — k)y^\knl{n — k — l)\ 



aoakaian-k-i^ 



4^3 E E 



6(n) 



xaQ^akaian-k^i + aian~k^iak)'h^ , 



(57) 
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where we have used e{k)/\k\ = 1/fc and changes of variables k n — k. In ( p7| ) 
we have obtained the two 0{h^) terms by commuting ak with ao and 1/A. At 
the next level the second term in (^5|) is calculated from [ao, a„_fc]2 as 



1 ■r-^ e(n) 3 



1 e(n) /I 1 1 

xaoakaittn-k-ifi^- (58) 
The third term in the parentheses in ( |55[ ) is calculated more easily as 

2^fc^n V\knin-k)\ 

^ 1 \ - \ - e(n) 2 

4^ ;^fer-. (" - kW\knl{n-k-l)\ 

1 e(n) /I 1 1 \ 

xaoaiUn-k-iak^^, (59) 

where we have used changes of variables k ^ n — k. 

Second, we calculate the second term in (|53|). As for the lowest term 
ao[[ao, an]2, ao]ih'^ , only the non-diagonal interactions give non-zero contribu- 
tions 



ao[[ao, an]nd2, ao]ndi?i'^ 

r e{n) / 1 1 1 

e(n) 1 /111 



e(n) V- /I 1 1\ 1 

«o I ao, — > , ( T + T + - ) :afca«-fc 

^' \j \nk\n — k)\ 



ndl 



X ao{akaian~k-i + aian^k-iak)^^ ■ (60) 



Other terms in the third term in ( |53| ) give higher-order terms. 
Collecting (|6l)-(^), we have 

ao[[ao, a„], ao] 

V V V — 

fc/o.n m=i ;#o,m 4A2m^|/cn(n - /c)?m(m - | 
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X ao{al^aiarn-i - ama]a^rn-i)^kan-k'h 



e{n) 



E 



1 1 



2^2 ,4:^ Vn ' k ' n-kJ y/\kn{n ~ k)\ 



1 



8^3 



k^O,n l^O.n-k 



1 3 2 

(n - k)yj\nkl{n - k - l)\^n k n-k 

3 



H ; 7 ao(afea;an-fc-( + aian-k-iak)fi 

n ^ k — i / 



(61) 



Appendix D 



The lowest-order term in [[[ao, a„], ao], ao] is [[[ao, a„]i, ao]i, ao]i?i^. In these 
commutation relations each [^,ao]i produces only higher-order terms. There 
are only three ways to get relevant terms. Choosing one diagonal interaction 
and two non-diagonal interactions for the three brackets, we have 

[[[ao, an]i, ao]i, ao]i?i^ 
1 e(n)e(n — k) 

'A3 2^ 2^ 



ni^o,n-k I" - k\^\nkl{n-k-l)\ 



4A3 

k=£0,n l=iO 
.11 1 1 

n I n — k — I 



+ T)^o{<iio-n-k-ia.k + akaian-k-i) 
k 



2 2 \ 

+ _ j^ o-oaiUn-k-iak + -^aoakaian-k~i)jfi^ ■ 

The other terms in [[[ao, a„], ao], ao] produce only higher-order terms. 



(62) 



Appendix E 



ao, CLnl, a 



"mj 5 CLm 



First, we calculate the lowest-order term in [[[ao, a„], am], aJn]' 



[[[ao, an]i, am]i, aL]i^^ = 

^ y- £(") rr 1 



rr 1 



j-aoa„, am 



2A 



nA 



(63) 



In the first term in (|6^) the non-diagonal interactions give higher-order terms. 
From the diagonal interactions we have 

1 



nA 



aoa„, am 



at 



19 



nmA^ 



2^2 



nmA^ 
1 



2^3 



(64) 



As for the second term in ( |63| ) only the non-diagonal interactions in [-j,a„ 
produce non-zero contributions 



e(n) 



rr 1 

2l 



e(ji)e(rn) 



my^lln^m — l)(n — m)\ 
n > 0. 



O-QdlO'm-lO'n-mTl' , 



(65) 



Second, we calculate the next-order term [[[qq, a„]i, a„i]j, a]„]fc7i'* (i+j+A; = 
4) . In these commutation relations only the diagonal interactions produce non- 
zero contributions. Paying attention to the commutation relations [-^,am], we 
have 

[[[ao, a„]2, a,n]i,alJi + [[[oq, a„]i, a„]2, aj„]i 



nA 



, a 



+ ) (1 + (5„,m)-^aoam?i ,n>0. 



(66) 



The rest of 0{K^) terms [[[ao, a„]i, Omji, a3„]2?i^ produces only higher-order con- 



tributions like {aoaaah ). From (p4D,(|65|) and ( p6| ) we have 

fi7 X! ^iii^o, an], am], al^] 
1 



6A TO 
m— 1 

C(2) .3 

E E 



12to2 ^|/n(TO - /)(n - to)!/!"* 



aoaiam-ian-mlT' 



1 



C(2) 



C(3) 
2nA4 



aoQnh ,n > 0. 



(67) 
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